Spectrum generating algebras for position-dependent 
mass oscillator Schrodinger equations 



C Quesne 

Physique Nucleaire Theorique et Physique Mathematique, Universite Libre de Bruxelles, 
Campus de la Plaine CP229, Boulevard du Triomphe, B-1050 Brussels, Belgium 

Abstract 

The interest of quadratic algebras for position-dependent mass Schrodinger equa- 
tions is highlighted by constructing spectrum generating algebras for a class of d- 
dimensional radial harmonic oscillators with d > 2 and a specific mass choice de- 
pending on some positive parameter a. Via some minor changes, the one-dimensional 
oscillator on the line with the same kind of mass is included in this class. The exis- 
tence of a single unitary irreducible representation belonging to the positive-discrete 
series type for d > 2 and of two of them for d = 1 is proved. The transition to the 
constant-mass limit a — > is studied and deformed su(l,l) generators are constructed. 
These operators are finally used to generate all the bound-state wavefunctions by an 
algebraic procedure. 
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1 Introduction 



During recent years, quantum mechanical systems with a position-dependent (effective) 
mass (PDM) have attracted a lot of attention and inspired intense research activites. They 
are indeed very useful in the study of many physical systems, such as electronic properties 
of semiconductors [1] and quantum dots [2], nuclei [3], quantum liquids [4], ^He clusters [5], 
metal clusters [6] , etc. 

Furthermore, the PDM presence in quantum mechanical problems may reflect some 
other unconventional effects, such as a deformation of the canonical commutation relations 
or a curvature of the underlying space [7] . It has also recently been signalled in the rapidly 
growing field of PT-symmetric [8, 9] (or pseudo-Hermitian [10] or else quasi-Hermitian [11]) 
quantum mechanics as occurring in the Hermitian Hamiltonian equivalent to some PT- 
symmetric systems at lowest order of perturbation theory [12, 13, 14]. 

Looking for exact solutions of the Schrodinger equation with a PDM has become an 
interesting research topic because such solutions may provide a conceptual understanding of 
some physical phenomena, as well as a testing ground for some approximation schemes. The 
generation of PDM and potential pairs leading to exactly solvable, quasi-exactly solvable 
or conditionally exactly solvable equations has been achieved by extending some methods 
known in the constant-mass case, such as point canonical transformations [15, 16, 17, 18, 
19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34], Lie algebraic methods [35, 36, 
37, 38, 39], and supersymmetric quantum mechanical techniques (or related intertwining 
operator methods) [7, 18, 19, 21, 22, 24, 34, 35, 40, 41, 42, 43, 44, 45, 46, 47, 48]. 

Another powerful tool used in standard quantum mechanics is that of nonlinear algebras, 
more specifically quadratic ones. For one-dimensional systems allowing exact solutions, 
such algebras may help us to understand the relation between the time evolution of classical 
dynamical variables and that of corresponding quantum operators, while providing a general 
method for constructing spectrum generating algebras [49] (see also [50]). In more than one 
dimension, they are a clue to classifying supcrintcgrable systems with integrals of motion 
quadratic in the momenta [51, 52, 53] and to solving the Schrodinger equation for such 
systems [54, 55, 56]. 
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In a PDM context, there has been no systematic use of quadratic algebras so far, al- 
though the presence of one of them has been signalled in a one-dimensional problem [43] . 
To start filling in this gap, we have recently considered the quadratic algebra generated 
by the integrals of motion of a two-dimensional superintegrable PDM system and shown 
how a deformed parafermionic oscillator realization of this algebra allows one to derive the 
bound-state energy spectrum [57]. 

In the present paper, we turn ourselves to another aspect of quadratic algebras, namely 
their occurrence as spectrum generating algebras, which we shall illustrate with the simplest 
example, corresponding to a harmonic oscillator potential. For a constant mass, it is well 
known (see, e.g., [58]) that all the states of such a potential with a given parity in one 
dimension or with a given angular momentum / in more than one dimension belong to 
a single unitary irreducible representation of an su(l,l) Lie algebra. The corresponding 
lowest-energy state is annihilated by the lowering generator, while the remaining states can 
be obtained from it by repeated applications of the raising generator. We plan to show that 
for a specific PDM choice, similar results apply except that su(l,l) gets deformed. We shall 
establish that a quadratic algebra approach provides us with a key to constructing such 
a deformed algebra, while allowing us at the same time to derive the bound-state energy 
spectrum. 

In section 2, we present the Schrodinger equation of a PDM d-dimensional radial har- 
monic oscillator {d > 2) and review its bound-state energy spectrum. The corresponding 
spectrum generating algebra is constructed in section 3. In section 4, we show how the 
general d-dimensional results can be applied to the one-dimensional oscillator on the full 
line with a similar PDM. Finally, section 5 contains the conclusion. 

2 Schrodinger equation of a PDM d-dimensional ra- 
dial harmonic oscillator 

Whenever both the PDM m(r) and the potential V{r) only depend on the radial variable r, 
the corresponding dimensional Schrodinger equation is separable in spherical coordinates. 
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On writing the radial wavefunction as r ^'^ ^^/^'0(r), so that the normahzation condition 

for ip{r) reads 

\^{r)\'^dr=l (2.1) 



+ Ves{r)]iP{r)^E'4;{r). (2.2) 







we end up with the radial equation 

did 



dr M(r) dr 

Here M(r) is the dimensionless form of the mass function m(r) = moM(r), we have taken 
units wherein h — 2mo = 1 and 

t7 / N T. / N {d-l)M' L(L + 1) 
Kff(r)^Kff(r)- ^ + 



2rM2 Mr2 

where a prime denotes derivative with respect to r, L is defined by L = / + (rf — 3)/2 in 
terms of the angular momentum quantum number I and Ves{r) is the effective potential 
that would arise in a cartesian coordinate approach to the problem (see equation (2.3) of 
[24]). 

Let us now consider a PDM d-dimensional harmonic oscillator, whose radial Schrodinger 
equation is obtained by replacing in the constant-mass one the radial momentum p. 



r 
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—id/dr by some deformed operator, tt^ = s/ f{oi; T)pr-\/ f{oi]r), where f{a;r) — 1 + ar 
and q; is a positive real constant. The result of this substitution reads 

(^^^^ + + ^) + l^V^ = E^"^i;^"\r) (2.3) 

which is equivalent to (2.2) with 

^^'^'^''^ = p{a;r) " (1 + ar^ 

and 

V;fr(r) = ^^^^ + l(a;2 - 8aV - « 

or 

Kff(r) = - Aa''[L{L + 1) + 2d]}r' - a[2L{L + 1) + 2d - 1]. 

Observe that the constant-mass limit corresponds to a — > 0, in which case equation (2.3) 
gives back the standard constant-mass equation. 



Sup CI symmetric quantum mechanical methods, combined with deformed shape invari- 
ancc, have shown [47] that the PDM Schrodinger equation (2.3) has an infinite number of 
bound states giving rise to a quadratic energy spectum 

El^l = a(^in^ + 4:n{L + l) + L + l + {'in + 2L + 3)^ n = 0,1,2,... (2.4) 

where A = |(« + A) and A = \/uj'^ + o?. In the same work, the lowest-energy wavefunction 
(for given V) has been obtained in the form 

^(^(r) = A/-Wr^+V-["+^^+')"l/^'") (2.5) 
where the normalization coefficient A/'g"^ can be easily determined from (2.1) as 



■2Q;^+§r(A + L + 2)' 



r(^ + i)r(^ + i) 

Some lengthy calculations along the same lines also yield 



^SW = ^Pi'"^''^^\t)4:i(r) (2.6) 



(---,i+-) 

where Pn ^ ^ (t) is a Jacobi polynomial [60] in the variable 

2 _ -1 + Q;r 
/ 1 + ar^ 



and 



/r(L+|)r(A + |)n! (A + 2n + L + l)r(A + n + L + l)^^/^ 



(2.8) 



M^i \ r(A + L + 2)r(A + n + l)r(n + L + |) 

Since in the constant-mass limit, the parameter A goes over to a;/2, it is clear 
that in such a limit the quadratic energy spectrum (2.4) becomes linear and given by 
-E'n.L — ^ (2?^ -|- L -|- I). Furthermore, the mere definition of e, combined with limit relations 
between orthogonal polynomials [60] also allows us to retrieve the results for constant-mass 
wavefunctions ■0n,L('"), depending on Laguerre polynomials [61].-*^ 



^Note that we obtain a phase factor (—1)" not present in equation (28.5) of [61]. This phase factor is 
consistent with positive matrix elements for the su(l,l) generators and with standard wavefunctions for 
the one-dimensional harmonic oscillator (see section 4). 



3 Spectrum generating algebra of the PDM d- 
dimensional radial harmonic oscillator 



In order to build a counterpart of the su(l,l) spectrum generating algebra obtained in the 
constant- mass case [58], it is useful to start from a quadratic algebra approach. It has indeed 
been suggested [49, 50, 54] that for a whole class of Hamiltonians, such as those for which the 
bound-state wavefunctions can be written as the lowest-energy one multiplied by increasing- 
degree polynomials in some variable t, there may exist an (in general nonlinear) algebra 
generating the spectrum, whose three generators are the Hamiltonian and the variable t, 
which are Hermitian operators, as well as their anti-Hermitian commutator. This algebra 
is characterized by a Casimir operator, which is some polynomial function of the three 
generators [49]. This is the approach to be followed in section 3.1. 



3.1 Quadratic algebra approach to the spectrum generating al- 
gebra 

Let us start from the Hamiltonian defined in equation (2.3), the variable t considered in 
(2.7) and their commutator, 

2 , L(L + 1) 1 



K + 



I 2 2 

+ -u; r 
4 



-Aia I 2-7rr + it 



(3.1) 



Prom the basic commutator [r, tt^] = i/(Q;; r), it is straightforward to derive the relations 



-f--l)+L(L + l)-l 



K. 



(3.2) 



16q;^ 



- (--1)-L(L + 1) 
a \ a 



showing that the operators k^\ k^^ and k'f^ generate a quadratic algebra. Its nature 
can be determined by comparing (3.2) with equation (3.2) of [49], defining the (general) 
Askey-Wilson algebra QAW(3) in terms of eight parameters R, Ai, A2, Ci, C2, D, Gi and 
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G2. Since in the present case, R = Ai = Ci = 0, we have to deal here with a quadratic 
Jacobi algebra QJ(3), characterized by the parameters 

"A /A 



A2 = -80; 
Go = -16q;2 



C2 = -16a 
A fx 



- 1 



a \ a 



a \a 
-1 1 - L(L+1) 



L(L + 1)-1 



(3.3) 



As — 4:A2Gi 7^ 0, this algebra is a nondegenerate one, i.e., an algebra that cannot be 
reduced to a Lie algebra by a change of basis. 

From equation (3.4) of [49], we get the corresponding Casimir operator in the form 

"A /A 



"A /A 



- 1 



a \a 



K. 



(«)2 



+ l&akf - 32q;2 



a \a 



- 1 -L(L+1) 



(3.4) 



K. 



Its eigenvalue can be obtained by inserting the explicit expressions (3.1) in (3.4) and is 
given by 



Q 



(a) 



16a' 



a V a 



L(L + 1) -2 



(3.5) 



Our aim now consists in constructing a positive-discrete series unitary irreducible rep- 
resentation of this algebra spanned by the Hamiltonian eigenfunctions ipl^iir), n — 0, 1, 
2, . . . , which will be a counterpart of the su(l,l) representation D'^ with ^ = | (-^ + |); 
obtained in the constant-mass case [58]. 

Prom the general theory developed in [49, 54], we know that in a basis ipp wherein the 
Hamiltonian, i.e., the generator k[°'\ is diagonal, the unitary irreducible representations of 
QJ(3) are given by 



Kf'ipp = gp+ittp+iipp+i - Qpttpipp-i 
where Ap, a^, hp and Qp are some real constants, which can be expressed in terms of the 
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defining parameters (3.3) and read 



Xp — a 



Ap{p + 1) - -(- - 1) - L(L + 1) + 1 

a a 



[16p2(2p-l)(2p+l)]-M2p 



A 



a 



+ L + 1 



2P---L 
a 



X 



2p + - - L - 1 
a 



2p 



a 



(3.6) 



^ = -|4p(p+ir'g-L-l) (A + L 

5(p = %ap. 

An infinite-dimensional representation of the positive-discrete series type is then 



characterized by the properties a^^ — and > Q \l p — p^ -\- n, n — 1, 2, Prom the 

exphcit value of given in (3.6), it is clear that, for generic values of A/ a and L, such 
conditions can be achieved in a single way, namely by assuming 



Po 



A 



a 



+ L 



(3.7) 



From (3.6) and (3.7), it results that the eigenvalues \pQ+n of k["^ in D^^ coincide with the 
energy eigenvalues (2.4), i.e., Xp^+n = E^^l, n = 0, 1, 2, . . . . 

Furthermore, if we reset Vpo+n ^ "0^1, ^po+n ^ bp^+n ^ and Qp^+n g^^l, 
the action of the generators K^""^ and X^"^ on the basis functions can be recast in the form 

-"■2 V^n.L — "n+l,L'A'nH 



where 



(a) 



-"■3 ^n,L — yn+l,L^n+l,L^n+l,L yn,L"'n,L^n-l,L 



'//(2/; + 2I + l) (^ + /; + !) (2^ + 2// 



(3.8) 



1/2 



+ 2n + L\ {^ + 2n + L-l){^ + 2n + L + l) 



(a) 



(A + 2n + L) (A + 2n + L + 2) 
A 



(3.9) 



= 4q! - + 2n + L 



a 



and T„ is a phase factor depending on the choice made for the relative phase of and 

^n-\ L- '^^^ fi^^^ equation in (3.8) can be reduced to the recursion relation for the Ja- 
(--- 

cobi polynomials Pn ^' ^ (t) and with the choice made in (2.8) for the normalization 
coefficients, we find that t„ = +1. 
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We conclude that the solutions of the PDM Schrodinger equation (2.3) can be derived 
by only using the quadratic algebra generated by the operators (3.1). To obtain from the 
latter the generators of a deformed su(l,l) spectrum generating algebra (and consequently 
a simpler construction of wavef unctions) , we shall need to build some ladder operators, 
generalizing the constant-mass ones. Before proceeding to such a derivation in section 3.3, 
it is worth considering the constant-mass limit of the quadratic algebra that we have just 
introduced. 

3.2 Constant-mass limit of the quadratic algebra 

Although appropriate for solving the Schrodinger equation (2.3), the basis 

(^k["\ K2"\ Rj^"^) of our quadratic algebra is not convenient to determine its a — limit 
because ^2"'* go^s over to the constant —1. To circumvent this difficulty, it is necessary to 
go from {ki''\ k^^\ k'f^) to a new basis {k'('\k^^\ R'f^). 
Let us set 

4 

^f) = l(l-A'("))-^(l + i^(-))=r^ 

i^f = ^{(1 - K^^Y^ ^3"^} = -2(2ir7r. + /). 
Observe that the inverse transformation reads 

^(a) ^^(a) ^ (i + Q,i^W)-i(_i + ai?(")) k'f^ =a{{l + aK^^'^Y\K'f^]. 

(3.10) 

Either from the commutation relations (3.2) of the first basis generators or by direct 
computation, we obtain for the second basis the commutation relations 

[K'f\K^^^] = 4{1 + - lea^A Q _ 1^ K^^\l + aK^^'^) 
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In the a ^ limit, it is obvious that these relations become hnear. It is then straightforward 
to show that the resulting operators Ki = lima^oKl"'\ i = 1,2,3, are some linear combi- 
nations of su(l,l) generators Kq, K^, K_, with commutation relations [Kq,K±] — ±K± 
and [K+,K^] = -2Ko. The results read Ki = 2uKo, K2 = {2/u){K+ + X_ + 2Ko) and 

Finally, on performing transformation (3.10) on the right-hand side of (3.4), the 
quadratic algebra Casimir operator yields, after some calculations, the relation 







\- 




a 





-1 +L(L + l)-2 



_ iSa^A (--1) + 8{i?i"\ + 12 (3.11) 



- 16L(L+ 1) 



+ 160«3^^") + 112a^i?(")'| (1 + aRi^'Y' 



Prom equation (3.5), it follows that the operator between curly brackets on the right- 
hand side of (3.11) vanishes. Since u)"^ = 4q;^^(^ — l), wc observe a close similar- 
ity between the first few terms making up this operator and the expression of the 
su(l,l) Casimir operator C = —K+K^ + Ko{Ko — 1) in terms of Ki, K2, K3, namely 
C ^ {Kl- ^uj'^Kl + 8{Xi, X2}) /64. We conclude that the substitution of a PDM for a 
constant mass has the effect of changing the constant C = \ {L + ^ [L — ^ into a function 
of r. 



(3.12) 



— [(2L + 3)(2L - 1) - lOar^ - 7q;V] 



3.3 Deformed su(l,l) spectrum generating algebra 

The purpose of this subsection is to construct a third basis {^kI^\ K^^\ K^^''^ of our 
quadratic algebra, satisfying the following three properties: 

(i) Kq"^ is proportional to the Hamiltonian of the problem, while K^^^ (resp. xi"^) is a 
raising (resp. lowering) ladder operator, which means that, up to some multiplicative 
factor, it transforms into V'i+i,L (resp. '4'^n-i,L) any n G M (resp. n G N"*") with 
the additional condition that i^i"'' annihilates iPq'l- 
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(ii) The operators Kq°'\ satisfy the same Hermiticity properties as Kq, K+, 

K_, i.e., ii^i")^ = ii-^") and K^^^^ = K^^\ 

(iu) In the a ^ hmit, they go over to the su(l,l) generators Kq, K^, K_. 

Prom the known action of k'^'' and ^3"^ on V'i"!) given in (3.8), we can construct some 
n-dependent ladder operators 

It is indeed easy to check that 

^+,n^n,L - (^n+l,L\9n,L + 5'n+l,L J ^n+l,L ^-,n^n,L " ^n,L[9n,L + 5'n+l,L j ^i-l.L • 

In (3.13), the quantum number n can be expressed in terms of E^l by inverting equation 
(2.4) and choosing the nonnegative root of the resulting quadratic equation. The result 



reads 



1 

n = - 
2 



- { ^ + L + 1] +Sn 



5n^\i^ + -(--l)+L{L + l). 
a a \a 



We can now eliminate the n dependence from by replacing E^j^ by the Hamiltonian 
H — . This leads to the operators 



^± — v-i- -t- -r -lu; 

where 



_ ^(a) _ 4q,x(")(1 if 5) + 4q;^^ ^ (3.14) 

1 ± () 



<5 = + A ("A _ 1 ) + L(L + 1). (3.15) 



Although such operators satisfy condition (i) referred to above, they do not fulfil the re- 
maining two conditions. 

We can get rid of this shortcoming by multiplying by some appropriate functions 
F^"^ [K^f^^ of the Hamiltonian. Since the latter are not univoquely determined by condi- 
tions (ii) and (iii), we may choose them in such a way that the action of k!^^^ on is 
the simplest possible. Let us therefore define 



^ ^l^At\b±\)\r-^ = ±^(5t 1)a/t^4"^ (3.16) 
^ 16A ^ ^ S 16A^ ^ ' \ St"^ 
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leading to the relations 



a 
A 



(n + 1) n + L + 



n+-+L+l n+-+- 
OL / \ a 2 



1 1/2 



r^{a) I (a) a 



n\n + L + -\ (n+- + L\ i n ^ 

a 



.1/2 



(3.17) 



2y \ a / \ a 

In (3.16), the factors ±-\/ (5 ± 2)/5 (alternatively ±-\/(5/(5 ^ 2)) are required by condition 
(ii) above, whereas the factors [5 ± 1) (alternatively (5 =F 1)) are optional ones having a 
simplifying effect on the matrix elements contained in (3.17). 
The definition of the third basis is finally completed by 

1 



(a) 



such that 



K, 



4A 



(a) 



-"■0 ^n,L — ^^^n,LVn,L- 



(3.18) 



In the a — > limit, equations (3.17) and (3.18) agree with the standard su(l,l) 
results K±'iljn,L{r) = [(n + | ± |) (n + L + 1 ± |)] ^''^ Vn±i,L(r) and -froV'n,L(r) = 

{En^L/2uj)^ri,L{'^), respectively. 

The three deformed su(l,l) generators k'^\ and satisfy the commutation 

relations 

[Af, At'] = ±^Al'>'(^ ± 1) = ±^(^ T DAf [At>, At'] = -"^ (2At' + ^) 

which can be easily checked by applying both sides on any ipl^^- Observe that for a ^ 0, we 
get a6/\ 1 and a/A — >• 0, so that standard su(l,l) commutation relations are retrieved, 
as it should be. 

The Casimir operator C^"^ of this deformed su(l,l) algebra can be written as C^°'^ — 
—K^^K^^ + /(Xg"^), where the function /(Xq"^) must be such that C^°^ commutes with 
and that /(i^^"^) ^ Ko{Ko - 1) for a ^ 0. The latter condition of course determines 
C(") only up to some constant term of order 0{a/X). After some rather lengthy calculations, 
we arrive at the result 

„2 



a 



8A2 
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leading to 



(3.19) 



Equation (3.19) should be contrasted with (3.12). 

In the appendix, it is shown how the ladder operators and can be used to 
fully determine the functions i/j^j^ in a much more direct way than those sketched above 
equation (2.6) and below equation (3.9). 

4 One-dimensional harmonic oscillator case 

The purpose of this section is to show how the results of section 3, valid for d > 2, can 
be extended to the one- dimensional harmonic oscillator on the full line. This implies, in 
particular, replacing the radial variable r (0 < r < oo) by x (— oo < x < oo). 

For a constant mass, it is well known that apart from the substitution r — > x, the 
Schrodinger equation for the standard one-dimensional harmonic oscillator can be deduced 
from the rf-dimensional radial one by setting either L = — 1 or L = 0. In the former 
(resp. latter) case, one gets the even-parity (resp. odd-parity) wavefunctions and corre- 
sponding eigenvalues, Vi/-i(?^)/a/2 i>2u{x), E^^^i (resp. V'i.,o(?^)/"\/2 'ijj2u+i{x), 
E^fi — > £'2,^+1), due to some relations between Laguerre and Hermite polynomials [60]. As 
a consequence, the single su(l,l) unitary irreducible representation D'^, = | (L-|- |), of 
the radial case gives rise to two such representations and -0^4 (with the same Casimir 
C = —3/16), for the one-dimensional case. 

The PDM Schrodinger equation 

TT = a/ f{a; x) p-y/ f{a; x) ^ ^ f {a; x) ^ 1 + ax'^ 



equivalent to 

d 1 d 



dx M{x) dx 
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admits a similar treatment exploiting the results obtained for equation (2.3), provided we 
distinguish again between the even- and odd-parity wavefunctions, given by 



'A 1 1 ' 



and 



'A 11' 



respectively. Here i/ = 0, 1, 2, . . . , t = 1 - (2//) = (-1 + ax'^)/{l + ax^), 



° "Ur(A + i)j 



_ /0Fr (A + 1) (A + 2.) r (A + .) ^ 



^ r(A + i)r(A + ^ + i)r(^ + i) 

^ / 2a3/^r(AH-2) \ A^li ^ /^ v/^r(A + l).! (A + 2. + l) T (A + . + l) \ 

^ lvir(A + i)j ^(.) 1^ 2r(A + 2)r(A + . + i)r(. + |) j 

and the corresponding eigenvalues are 

in both cases n — 2v and n = 2z/ -|- 1. 

There exists a quadratic spectrum generating algebra, for which we can construct three 
sets of generators X^^^), {K^^\ K^^\ R'f'^) and {K^°'\ K^^\ K^^^) , analogous 

to those built in section 3. The only differences lie in the substitutions r ^ x, Hr t^, 
L[L + 1) — > 0, and in the very important fact that there are now two distinct unitary 
irreducible representations instead of a single one. This can be seen from the counterpart 

al = [16p2(2p - l)(2p + 1)]-^ {2p - ^) - ^ + 1^ {2p + (^2p + ^ - 1 

of the similar quantity defined in (3.6). The conditions a^^ = and > if p = po + 

= 1, 2, . . . , characterizing positive-discrete series representations are indeed satisfied 
now by two distinct values of po, Po = | — l) and Po = corresponding to L = — 1 and 
L = in (3.7) and to which we can associate Xp^+v — and Xp^+v — -E'2°+i) respectively. 

Since, after these observations, it is straightforward to transpose the results of section 
3 to the one-dimensional case, we are not going to detail them here. We would only like to 
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mention that the action of the deformed su(l,l) generators on the wavefunctions reads 



i^i^Vi"^^) = ^ - 1) + 2- - 2^ + 2- - 1^ 
4A[ \ a J \ a ) 



leading to the standard su(l,l) results Kot/jn{x) = | + ^)'4'n{x), K^ipnix) = |[(n ± 
l)(n + 1 ± 1)]^/Vn±2(a;) in the a ^ limit. 

5 Conclusion 

In this paper, we have highlighted the interest of quadratic algebras for PDM Schrodinger 
equations by constructing spectrum generating algebras for a class of d-dimensional radial 
harmonic oscillators with d > 2 and a specific PDM choice, depending on some positive 
parameter a. We have also shown how minor changes enable the one-dimensional oscillator 
on the line with the same type of mass to be included in such a class. 

For these quadratic algebras, we have considered three different sets of generators. The 
first one {k["\ ki'^l K^^""^) has allowed us to prove the existence of a single unitary irre- 
ducible representation belonging to the positive-discrete series type for c? > 2 and of two of 
them for d = 1, as well as to obtain the bound-state quadratic energy spectrum. 

The second set [k[°'\ K2°'\ K^""^) has provided us with an exphcit demonstration that 
the quadratic algebra considered here gives rise to the well-known su(l,l) Lie algebra gen- 
erating the oscillator spectrum in the constant-mass limit, i.e., for a — * 0. 

This correspondence has been studied further by constructing a third set of operators 
{Kq"'\ -f^i"^), which go over to the standard su(l,l) generators {Kq, K+, K_) for q; ^ 
and may therefore be termed deformed su(l,l) generators. All the bound-state wavefunc- 
tions have finally been built by using the lowering and raising generators, xi"^ and 
respectively. 

Some interesting open problems for future work are the extensions of the present study 
to other exactly solvable PDM Schrodinger equations either with the same potential but a 
different mass or with both different potential and mass. 
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Appendix 

The purpose of this appendix is to prove equations (2.5)-(2.8) by using the deformed su(l,l) 
algebra introduced in section 3.3. 

Let us start with ip^^^lir), which, according to the second relation in (3.17), is annihilated 
by K^^ or, equivalently, by A^^\ Equations (3.14) and (3.15), together with (3.1), yield 
the first-order differential equation 



whose solution can be written in the form (2.5). 

The excited-state wavefunctions ■^I'lI'")' = 1, 2, . . . , can now be determined recur- 
sively from '0o'l(^) employing the first relation in (3.17). When combined with definition 
(3.16), the latter yields 



^Sl,L(r) 



( 2n + - + L + 2 
16q; \ a 



X 



n+l)in + L + 



1/2 

2n + - + L + 3 " 
a 

n+ - + L + 1] (n + 
a 



f ^ 




n + - H 




\ a 





1-1/2 



(A. I) 



X (2n+^ + L + lj 4"VS(r). 



Let us now make the ansatz 



/vr(«) 

^SW = -^V^S2(r)Pn(^) 



(A.2) 



where P„(i) is some nth-degree polynomial in t, such that Po(^) — 1- On inserting (A.2) in 
^+ Vi"l(^) ^-iid using equations (3.1) and (3.14), we get 



4"Vil(r) 



-8a 



1 1 2n + - + L + 2 ) (l-t^)^ 



— n 



'0,L 

A 
a 



L + 1 



X 



a 



-L-1+ (2n 



A 



L + 2]t 



Pn{t) 



which, according to (A.l) and (A.2), should be proportional to ij/^°'l{r)Pn+i{t). This clearly 
identifies Pn{t) as the Jacobi polynomial Pi^''^\t) with P— ^ — ^, ^ — L+^, because the 
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latter satisfies the relation 



|(2n + /3 + 7 + 2)(1 - t^)! _ + /3 + ^ + _ ^ + (2n + /3 + 7 + 2)i]| 
X piM(t) = -2{n + l)(n + /3 + 7 + W 



(A.3) 



obtained by efiminating P!^i\t) between the Jacobi recursion and differential relations (see 
equations (22.7.1) and (22.8.1) of [60]). Hence equation (2.6) is proved. 

Finally, on combining equations (A.1)-(A.3), we arrive at a recursion relation for the 
normalization coefficient 

•A/"i+i,L / (n + l)(n + A + L + l)(2n + A + L + 3) ' 
M^l l^(n + L+|)(n+A + l)(2n+A + L + l) 

whose solution is given by (2.8). This completes the determination of the wavef unctions 
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